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ABSTRACT 


The Coppens-Sanders perturbation solution for the one-dimensional 
non-linear acoustic wave equation with dissipative term Geecounire 
the viscous and thermal energy losses encountered in a rigid walled, 
closed tube with large length-to-diameter ratio was extended to include 
sixth order terms. The solution was then investigated to determine 
the region of validity. Computer programs were written to evaluate 
and graph the resulting waveforms. Available experimental results 

ee 

were compared with the theoretical predictions and good correlation 
was found to exist in the region of low Mach numbers. This agreement 
was found to gradually deteriorate as the Mach number was increased. 


A Fourier synthesis approach is also presented and the leading terms 


of the first ten harmonics are derived. 
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operator for body forces (see Eq. 9) 
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1. INTRODUCTION 


The thermodynamic equation of state of a fluid for reversible 


acoustic — can be written as 
arp) hams sae (/ -p.)* (1.1) 


or; for an ideal gas, 


‘wit he ty \" 
mth, 





Gleg2) 


where J is the ratio of specific heats. The continuity equation for 


one-dimensional wave propagation in Lagrangian coordinates is 


pir ee (aes) 


ed 


The force equation in the absence of dissipative mechanisms has the 


dt __ of. 


J A ps (1.4) 


form 


It is iS sali by straightforward combination to obtain the equation 


a sd 7 g) fi 
Of /), MF "OF _ (1.5) 
D5. (152 7° 
for a fluid obeying Eq. 1.1 or 


2 - 
ot. ae (14 e! oF | =a) (1.6) 


dae i 


for an ideal gas. 
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If interest is restricted to acoustic processes foreuaien oy <1, 


Eqs. 1.5 and 1.6 may be approximated by 


=~ (2 # . idea 


which (upon differentiation with respect to time) becomes the non- 


dissipative wave sy? 


cr lu - fs z)° (1.8) 


p> tel Gy) om SS Py 


ll = die «oe 
aa” ae at? 


where 


A more general form of Eq. 1.4 is 


gu __ a 
fo ve = ape a (1.10 


Mier is an operator which generates all forces, other than that 
arising from the gradient of the pressure field, which are active in 
the system of interest [4]. If Eq. 1.10"1s Wsed instead of aq. 1.4. 


Eq. 1.7 becomes 


=) rad + (Gz a “<2 


In the Navier-Stokes case 


eo 
4 = (0+ ej a 


where % and %, are the shear and bulk viscosity coefficients, respectively. 
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An important contribution to the solution of Eq. 1.11 for the 
case of Eq. 1.12 was made by Fay [5]. He investigated the changes 
in the steepness of the wavefront in a periodic, finite-amplitude 
plane wave of infinite extent in a viscous medium. By means of 
Fourier analysis techniques he was able to show that the nonlinearity 
of the pressure density relationship results in accumulating distortion 
of the waveform with shifts of energy from the lower to the higher 
frequency components of the wave. Because the classical Navier-Stokes 
absorption coefficient is proportional to the square of the frequency, 
the higher frequency components are attenuated more rapidly. Fay's 
solution is valid in a region far from the source where the rates of 
harmonic growth due to nonlinearity and the rates of attenuation due 
to viscosity are tending to balance leading to a waveform of relatively 
stable shape over distances corresponding to many wavelengths. 

Fox and Wallace [1] investigated the same problem as Fay using a 
graphical analysis technique. Their result is essentially the same 
as Fay's. Keck and Beyer [6] developed a perturbation analysis technique 
for periodic plane progressive waves of infinite extent in a viscous 
medium with which they generate terms through sixth order. The wave 
equation is the same as that treated by Fay, but the solution is valid 
only near the source. 

Weston [20] has presented a linearized wave equation for the 
propagation of monofrequency sound in tubes. He assumes that the wave- 
fronts are planar except near the walls. The viscous and heat-conduction 
losses in the boundary layer are assumed to be the dominant dissipation 
mechanisms in the fluid. Except in this boundary layer, the wavefronts 


are independent of a radial coordinate. This suggests that the boundary 
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layer loss can be replaced by an equivalent absorptive process active 
throughout the volume of the cavity [4]. More will be said about this 
later. 

Saenger and Hudson [14] present a simple description of periodic 
shocks at resonance using a two part solution. Their solution is 
based on the assumption of the applicability of the linearized acoustic 
equations everywhere except in the region of the shock where the Rankine- 
Hugoniot shock conditions are assumed. They consider the case in which 
the piston oscillates at the fundamental frequency. Compressive viscosity 
is ignored but the effects of shear viscosity and heat conduction in the 
boundary layer are considered. Their solution remains finite only because 
of these considerations. 

Betchov [13] also postulates the existence of the shock wave and 
constructs a solution at resonance based on a continuous and a discontinuous 
part. For an inviscid fluid the amplitude at resonance is found to be 
finite and to be determined by nonlinear effects. The effect of wall 
friction is discussed and it is suggested that this could modify the 
solution significantly. 

Weiss [7] has applied finite-difference techniques to the nonlinear 
inviscid acoustic equations and has shown that repeated reflections at 
the rigid boundaries tend to promote the development of a discontinuity 
in the velocity profile. The absence of any dissipative mechanisms, 
however, eliminates the possibility of any steady-state standing wave 
patterns. 

While most investigations have been confined to traveling waves, 


a few have considered standing waves. 


Me 


Solutions of the nonlinear wave equation have been obtained for 
the case of finite-amplitude standing waves in a closed tube, where 
the tube is driven at one end by a piston source. This work was done 
by Keller [21] using the Lagrangian formulation and assigning an 
unrealistic value of -l to the adiabatic exponent ¥ . Due to neglect 
of any dissipative mechanisms, his expression for the particle velocity 
becomes infinite if the piston frequency equals a natural frequency 
of the tube. 

A detailed theoretical analysis based on the nonlinear acoustic 
equations including the effects of compressive viscosity and shear 
viscosity in the boundary layer has been done by Chester [8]. He 
successfully predicted asymmetries in the finite-amplitude standing 
waves. Unfortunately his results are not easily compared to experi- 
mental results. 

Coppens and Sanders [4] developed an extension of the Keck-Beyer 
perturbation approach wherein wall losses as well as bulk losses are 
included. Their solution is applied to the case of finite-amplitudes 
standing waves in rigid walled cavities. The results yield information 
concerning the amplitudes and phases of the Fourier components of the 
waveform, and indicate the importance of the type of absorptive process 
on the resulting waveform. 

The purpose of this thesis is to extend the Coppens-Sanders 
perturbation approach. In part this is accomplished by a direct 
application of their method to obtain sixth-order terms. Information 
concerning the response of the cylindrical cavity to the sinusoidal 


excitation is extracted from this solution by some algebraic manipulations. 
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A Fourter synthesis approach to the problem is presented and from it 
the leading terms of the first ten harmonics are derived. Information 
similar to that obtained from the perturbation approach is acquired 


and both results are compared to the experimental results obtained 


by Beech [19]. 
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2. THEORETICAL DEVELOPMENT 


In this section both a perturbation approach and a Fourier-synthesis 
approach to the solution of the one-dimensional non-linear, dissipative 
acoustic wave equation are presented. These approaches are formulated 
for a rigid-walled, closed, cylindrical cavity with large length-to- 
diameter ratio. The solutions obtained are for finite-amplitude standing 
waves. Manipulations are also presented which enable each approach to 
yield resonance response information which is convenient for comparison 


with experimental results. 


Perturbation Approach 


Assume a perturbation series ok ee such that the n-th 
nn 
term is of order n-l in the Mach number. Substituting this into Eq. 1.8 


and collecting terms of equal order yields the set of equations 


(= : D) uy = D, oo 





where 
Ey 
p d 77 oF’ Bee Wy 
io ) 
D, = me, (2.2) 
O n< | 
Tf Uyn= 2s nj , where the summation over j is understood to 
d 


encompass only those frequencies included in each Dy , Eq. 2.1 becomes 


2 (OL a Dj ) un; = Dy (2.3) 
7 


=o 


In order to utilize this formulation, it is first necessary 
to obtain the form of Dj. For the case of finite amplitude sound 
in a duct, Weston [20] has presented a wave equation valid for the 
propagation of monofrequency sound, Weston's equation may be written 


in Eulerian form as 


oY 
ee + § (i-!) = = O gues 
where 
-/o| |! ~| Pe! 
§=(S)(20) “[v4+C-F )V” 0.5) 
and 


It 


(uniform) cross-sectional area 


perimeter of area S 


(angular) frequency of the propagating waveform 


YY © maee 
il 


kinematic shear viscosity of the fluid 
). thermometric conductivity of the walls 
is Ya 
This may be rewritten as 
3, /f) ot - s iC (2.6) 
Lal pe (5) at 7 we | 
gx 
When dissipation processes are included, it is usually assumed, explicitly 
or tacitly, that they are weak. Under this assumption, and for small 
Mach number, the derivatives (Wa yt and (Say )t may be interchanged. 
This also holds for the pair (Sp a and (Yt)x. Thus Eq. 2.3 must be 


solved with 


3 2 
Ys — TP abe 
Dj = (i) s808 ( da™ (2.7) 
where di represents the value 6Ee tes a frequency a; : 
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The final result is the approximate wave equation 
3 eZ 
ci* + ie Sato 
ites C7 OS aFn-é 
Pees =— 


Eq. 2.8 was solved [4] for a uniform tube of length L excited 


(26) 


at end a = 0 by a piston driven with constant acceleration A, at a 
frequency @ , and terminated at a = L by a rigid cap. In this 


solution, the boundary conditions were given by 


uy, (o,t) = (AYa) tor (wt - 6) 
(2.9) 
Un Ge = Uy, (9t) = ding (t,t) mance wacy) > | 
where o will be defined later. 


The lowest order solution of Eq. 2.8 is 


i = Ce (4) (2.10) 


where 


wih = Pe: -Q 
Uy, = Ao o( (et 4) € Ea) e ig ea) 
- ee) 


(2ecuhal:) 
hole =m 
e ~ sae, 
and k = k -i&® is the complex propagation constant. 
The solution for Wi may be rewritten as 
Uy, = UI, b[eosk (4-2) sin h(L-a)] cca wt 
(2.12) 


+ tinh x(L-a.) CE2- kék-a) | ti. at } 


where 


= (2) (ah wl + whl)? (2.13) 


an 


and 





rai el 
/ a Taw — (2.14) 
Tuk CL 
Eqs. 2.12 and 2.13 will be rewritten in terms of Oy , the 
resonant frequency which maximizes the pressure amplitude Pui ' 
associated with U1; at the rigid end of the cavity. 
Z 
yp ee 
—_a = 7 
C jis ma (2.15) 


ww = wp + AW 


The integer m represents the normal mode most strongly excited by the 


GC: 
input frequency. For ol <« | and AG <7Uf and O, are 





77 ZA 3 
= EN 2.16 
HW mr, con ( ) 


G: = Tons’ (~ fees ) GAP 
Jj op 


From this point on, higher-order solutions are obtained by first 





making the approximation 


Loo 
Ay = Ly K (L-a) ecu wt (2.18) 


which is Eq. 2.12 with terms of Wh omitted, The procedure (presented 
in Appendix A), is then merely an iteration process. 

In view of later discussions it is advantageous to point out that 
Eq. 2.18 may be rewritten in terms of the associated acoustic pressure 


by merely making the substitutions 


“L— TY po0o UU — fat 


dw > Oct LeL— Mw 


C279) 


ZZ 


This is equivalent to ignoring terms of order Yp in Eq. 1.10, and 
converts the calculated variable from velocity to pressure, a more 
commonly measured quantity. 

In terms of the pressure, the solution presented by Coppens and 


Sanders is 


Yip, = coe Ke} aivot -H(B2, a. 
+4 (SE )*/4,3.2,3),-£ (42-23), 
t(U4 22) +z (42-2) Je f 
toor 2K (L-a) | $ (Bia), - UAT e24, 
+(422h] tof 
+ eo FR L-a) ve (EVs fF) fe23%, 


24) 


i (2,334); i (OAS, OF 


yee a: 


+2 (22,3), [4+ } 


teow # Klbra WA (FY [i254 lytg 284) ]-f 
a5 ca SK «) 4 (BY (2 34 5), 
+ pees *3(2255)s/4.. ; 


pees 


23 


where 


(4,b)¢ - /? 7 yon (gut 46,744) (2.21) 
- 53 
Hy = [Cby-& 7 2aypy)e +h? f° ae 


ie at a y, 
On = ~ Cot [OV Syd + ary, fy ney (2.23) 


On /o (2.24) 


and M is the peak Mach number of the first order solution. 
Using the same scheme, the sixth order terms were derived. 


The solution for Pe which is to be added to Eq. 2.20 is 


ie = Con 2K(L-a)! 2 (BY faz 0 2 3 ¢/y (27-259 


(22234), + (22233), 


24 


oT - 
- 002 4 Kea), HF) /@. 44,5), 
AZ (23445)4 +3 (223,45)4 
r(23344)4 +d (223 44)s 
4(22 2349) 4(12234)y 
#3(4225,4)4 +4(22234), 
/ 
Seta is (1323 ¢h]+~-| 
5 
t Cor bK(L-a } HA Soe 
Ay g 2245 6), 2. Zz 235 by 
/ 
*2 (24234, +3 (2z24h)i 


£02254] 
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The entire perturbation solution herein developed, which is the 
sum of Eqs. 2.20 and 2.25, will be discussed in some detail after 


the Fourier synthesis approach has been presented. 


Fourier Synthesis 


Beginning with the approximate wave equation Eq. 2.8 and assuming 
the total solution is composed of the known classical solution §% 


satisfying 








bf 2 
yer+s ast _ 5; ey : 
< 4] da&ot da 


and some other portion ¢° defined by 


=O (2.26) 


sh? (2.27) 
then 
CC eee G of ay aie © 
| godt | ga™ 
. fie Or< , oe 
Z Oa Ja” Ja 





(2.28) 


Define the phase velocity Ci to be G7 =0.2 (Hd) and Eq. 2.28 may 


be written 


S ca 
(zi = ‘e Gt (A) 25 at 


y+ a os ae )* 
= "7 vce ag 


C229) 


Assume 5 may be written as a Fourier series 


pees 
S es Coa NK(L-a.) cot (nwt + Pn ) Q2E30) 
=O 
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and 


ok =X Cor K(l-«) Cor me CFET®) 
my 


and let 


a “ 57 Aloe) co+(nust ip [> ) (2.32) 


where An, and iV. will be obtained later. Combination of Eqs. 2.29 


through 2.32 yields 
= E i 
L(z/- (2. *[% bi WK(L-%) Cor brut t Da 4) 
=n dow MCLEE) divs (ot + bn?) (2.33) 


2 tn (2) cor (nwt +/m ) 


We impose the conditions that 
Ay (&) = iowcee nk (La. ) (2.34) 


and 


Anta) < nK sin nKE-2) 


so that Eq. 2.33 becomes 


(2 >) oe) | cbsatt bn) gn‘) - 5 ssfuat 1 po’) f Xo 
= = nen 7 as) 


Introduction of the phase angles 0, of Eq. 2.17 allows Eq. 2.35 to be 


(2225) 


written 





Ye Ain rat + por - Oy) = a ae Bn corfnut +i) (2.36) 
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In order to simplify the manipulations to obtain Be and /;, 


define 


X) Cor K (L-c) Corlut + d,) = X Coa KtL- 2) Cos(wt # g, °) 


+X Cot K(L-@ ) Cone 


Xn Cor nk(L-«) Corfiwt t Py ) 


=X, Cor nK(l-a) Cot brute da?) , Nn ae 


Thus 


os 
= 2, Xn Cor MK(L-a) Con (nt + Pn ) 
n= 


2B, co ‘eal +i) 
=), XG Xone cor (nwt + Bay 5 - 4; ) 
-/ 


t 2% Xne4 cod (nut + Py-4 +d; ) 


Combining Eqs. 2.40 and 2.36 gives 


_4 da 


Yr on On, ar Lin l(nat4 g,, - On ) 


ms gi Co2 (nut Pant - 47) 
ae 


‘% 2G Xo Coe hut +$n3 +4) 


ey 


28 


(2 


(2% 


(2; 


CZ 


(o7) 


~41) 


It is convenient to shift to complex notation at this time. Let 


X, Cor F, a ee Cot 7 + Xe 





(272) 
x hin. Bi aa Ai yu. @ 
so that 
ra) cg,° 
ae ae (eet3) 
and 
2 5 page 
ie s \,. : on > (2.44) 
Then Eq. 2.41 — 
« i (ti’- 6 -% 
-+= a ms ) 
(2.45) 


ty; x (45-3) 
= Kt € 
LS ie 
= 
If the equivalent pressure , aw is desired, make the 


32 Yh oi (Gr; +65) 


substitution 


ei = - Yo,00 (2.46) 


where M is the amplitude of the Mach number for the particular harmonic. 


Define 


pL, on 


so that Eq. 2.52 becomes 
oO = 
Mey oi Con’ om = es) 


Dy 
<7 ; =e 
= Zo! Mai; € (Putz - $7) — 
d- i =i 
+2 2. MG Mang @ OI j +97) 
ig { 
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Examination of Eq. 2.48 and Eqs. 2.20 and 2.25 shows that the two 


methods are equivalent since expansion of the term 
n-1 


Z ZL Mz IT 4-5 e Fm 3 * $7) 


yields the leading vem of each harmonic derived in the perturbation 
scheme, while the other portion of Eq. 2.48 contains the correction 
Eerms . 

Recall that in both the perturbation approach and the Fourier 
synthesis approach the leading term of the fundamental was calculated 
to be li pa, K(L-a) Co. wt - This is of order one in the Mach number: 
It is also a known experimental fact that the fundamental remains 
essentially unchanged with respect to Mach number. This suggests that 
the phase angles associated with the corrective terms might cause the 
corrective terms to sum to approximately zero. It then becomes 
conceivable that a solution composed of only leading terms could 
yield significant results. In order to investigate this possibility, 


Eq. 2.48 was used to obtain the leading terms of the first ten harmonics. 


The result is 


% = Z Cor Ka) | 2. wt} (2.49) 
+ 5 Cow 2x(kea)t (SP) a, 
+3 cor 3K(L-0) {(Sf)%2)s§ 


pte ence) (QE) fess 


30 


+ £02 4/4] 
+ cor 5K (La) { (F) /@3 457s 
te (2245), Slezseie]} 
+t orb Kha) (HY fess y, 
+f (2245% +4 (224556) 


4522546) +f (22246), 


Pa 


aN 


ae nS, oJ} 
13 Con 7a) (SF) 234567) 


/ _ = 


7 


+z (823461), + 2 (222461, 


5 ae 


uv e | 
tg (< ee 67), aS. (22 ane? 


My 


/ / 
+3 (222957), +4 (222357) 


oS ee: 


i f / 
4 3 (2295 47), + B tess 


oe. 


m 

+ 3 Co4 sKcLal{(F4) [(2145:624) 
i - 

#3 (229567 Dp +3 (2235678), 

"222246199 + $ (2224 7, Oe 


+ $22 $2678), +5 (2239578) 


yA vi } pe? 


3 (2224579), a d (2235, 72) 


(2227 


2 C2 55474s + Z (2223478) 


/ a), 1 
oy ee £56 C/o 
! : / 
4g (FES 556E) + £ (222 3468)q 


/ 
76 (2222.46 8), * 3(22237 68), 
ioe 58 tee 
yy Sahay , }, & (222545 Je 
y a 
+g (2225358) + 5 (22 33-948) 
+ $G225498), 44 (222.446) y, 
B 
#7 Cee IEC éa), ( of) i Z54507EN, 


L 
tg (2245078) + M2235 76%), 


#3 (2546789), + £02246 769), 


2G: 


45(22346781)q + $ (22345781g 


/ Pa 
+3 (2234562 T]y + 3 (222456 PEs 
+g 22255689), +f (ZA I46 S44 


/ 1 
(22224689), 4 3 (222 33609} 


oe) De Fe 
f - 
4 t(22F7S4¢5EF al 2 
errr 4 * 3 (433 749E% 


+— Z (22253589) + £ (2233 334481), 


| ee dae A227 7 


Saf fj} 


#4 jf (222556) 77), + ‘ Z (2223467 9p 


/ 

+ j¢ (2222457 I) 4 # 2 (2225379), 
iL 

Ty (22254 574), + be + (22224579 )g 


pyeg s 


i 
43 (22225579), + 4 (22277479); 


23 


/ 
+4 (2222i479)9 + 3(2353456 9 
t (2223456 I) + £ (22833559) 


ae: 


/ oft 
t 7 (22253 469y * 7 (3332346 9)o 


od. 
+3 (Z225336 Ip + 322334959 Jy 
/ 
4 B(2225 4450), + F (2223345 9) 


t 3 (222344 55 + ; 3 (23324459), 


7 
+3 Cor (0K (L-a} iF ) B34), 
/ 57 
+ 7 (224 56/8 7/6}, + 2 £ (2233 97 EG 16) 
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where (a,b)g has been given in Eq. 2.21. 
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Cavity Resonance Response 


Investigation of the behavior near resonance of each approach 
leads to the desirability of obtaining the response of the cavity 
for each harmonic as the fundamental is swept through resonance. 
ft is obvious from the form cr Bas-)2.2/, 2.32,1sand 2.55 thareborn 
approaches yield results which may be written as sums of terms 
having like factors of K(L-a) andwWt. It is also evident that each 
of these harmonics is heavily dependent upon M, the Mach number, and 
Aw. 

The algebraic process necessary to bring each of the two solutions 
to a form which conveniently yields the desired information will be 
derived in a general form and then applied to the perturbation 
solution. It is then easily extended to the Fourier synthesis solution. 


We begin with the relationship 


yw Lin (fut + $2) = A to Gat + ). (2.50) 


| 

At t = 0 Eq. 2.50 becomes 

— | a (2.51) 
We Qo hw di - A Lin g. 

Cat 


The derivative of Eq. 2.50 evaluated at t = 0 yields 


7) 
UY ¢ a 
Le Le Cor Oe = Noor p (2.52) 
Ca 
which when combined with Eq. 2.51, yields 
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If Eq. 2.51 is divided by Eq. 2.52 the result is 
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The sum of Eqs. 2.20 and 2.25 may be written as 
nn 


Pip, = pcos jkoolZ. dr die yat rg) 


(2255) 


4 A ove ghee) di. jute g. 


for each harmonic, where j indicates the harmonic number and n 
indicates the number of terms in that harmonic. It is evident 
that the amplitude of each harmonic may be obtained by applying 
ge 2.03. 

In applying Eq. 2.53 to the perturbation approach the magnitude 
of each harmonic was normalized to approximately one by the maximum 
magnitude of a factor chosen from the leading term of that harmonic. 
This normalization was done to yield a set of curves whose maxima are 
of order unity. For the fundamental this normalization factor is one. 
The method of choosing the normalization factor for higher harmonics 
will be demonstrated for the second harmonic. The leading term of 


the second harmonic is 


(Fhe, £4 of htw (2ut +62), (2.56) 
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and the leading term of the second harmonic may be written 
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Since 


Ay “= em G,, (2.60) 


Eq. 2.59 may be written 


ofp Caan Caz G2 Luft 42 ) (2.61) 


2m Tl p Co dy dz 





H has as its maximum 


#5 Irae — a (2762) 


and the normalization factor for the second harmonic is then chosen as 


Ae oe (2.63) 


Lm Kt, Co 
It should be noted that these normalization factors are not the maximum 
amplitudes of the various harmonics, neither are they the maximum 
amplitudes of the leading terms of each harmonic; rather, they are 
factors chosen to yield maxima of order unity. 
These factors are given the symbol MAXA where n refers to the n-th 


harmonic so that 


MAXA., = Ap + (2.64) 
4 Lm A by Co ds d2 


By combining Eqs. 2.53 and 2.64 it is possible to define o. as the 
normalized response of the n-th harmonic. It should be noted that A 
defined in Eq. 2.53 is a function of the Mach number and AW , 


Thus 


Q), (M Aw) = AY OKA + he 


where n is the harmonic number. 
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3. APPLICATIONS 


Associated with each of the solutions formulated in the previous 
section are three computer programs. These programs are concerned 
with the pressure waveform, the Q-curves, and the phase information. 

The perturbation solution has the programs FINAMPI, QCURVES, and PHAMP 
connected with it while FINAMPIZ, QCURC and PHAMPC are related to the 
Fourier synthesis. 

The first attempted utilization of the perturbation solution 
consisted of writing a computer program involving all terms through 
sixth order. This program, entitled FINAMPI, was designed to compute 
and graph both velocity and pressure waveforms. FINAMPI was first used 
to graph the waveforms contained in Ref. [4] and has been used since 
in determining for what Mach numbers the theory begins to deviate 
Significantly from the experimental results. The graphical output 
from the computer is contained in Figs. C.1l through C.4. 

The program designed to compute and graph the cavity resonance 
response for the perturbation Polueson is entitled QCURVES. This program 
is a direct application of Eq. 2.65 and has been used to obtain the 
curves contained in Figs. C.5 through C.22. 

PHAMP, the third program associated with the perturbation solution, 
was designed to compute the phase angles of each harmonic as the 
fundamental was swept through resonance. It is a direct application 
of Eq. 2.54 and resulted in the data contained in Table 4.2. 

In an attempt both to determine the importance of the corrective 
factors to each harmonic and to test the validity of the Fourier synthesis, 
computer programs involving only the leading terms of the various harmonics 


were written. Unfortunately a program involving all ten harmonics could 
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not be brought to fruition because of time limitations. A second 
program, FINAMPIZ, which uses the leading terms of the first six 
harmonics has been used to graph the pressure waveforms also contained 
ingkigs..C.1 throuwehgG. 73 

Eqs. 2.65 and 2.54 were also applied to the Fourier synthesis 
yielding the programs QCURC and PHAMPC. QCURC computes and graphs 
the Q-curves for each of the harmonics while PHAMPC computes and graphs 
the phase dependence of the individual harmonics on AW, 

It should be noted that of the two programs concerned with phase 
dependence, only PHAMPC graphs its output. This was done because the 
Q-curve programs indicated the Fourier synthesis approach produced 
curves which were in much better agreement with the experimental results 
than did the perturbation approach. The graphical output of PHAMPC is 
contaamed in Fig. C.23. 

In this section only a superficial description of each program 
was given. A detailed description of all programs is contained in 
Appendix B, and a more complete discussion of results follows in 


Section 4 


42 


4. RESULTS AND COMPARISONS 


In this section the theoretical predictions obtained from the 
applications described in Section 3 will be compared to the experimental 
results obtained by Beech [19]. 

Before proceeding to a discussion of these results it is necessary 
to attempt some qualitative evaluation of the apparent properties of 
the perturbation solution. The amplitudes of the (n + 1)th order terms 


tend to behave as [4]. 
n 
n awe 
27M J | es! 
jr< 7 
It is clear that this expression can diverge for increasing n. This 
divergence, which occurs because of the dependence ine rm , would 
mot occur if bed since the left hand side of Eq. 4.1 would then 
y 
be bounded by (AS) and all series would converge for “ps 68 
: ee : iM Fs 
The Mach number which satisfies the expression Wade / and would 
therefore appear to be related to shock formation is M = 0.02. This is 
within the same order of magnitude as the experimentally measured Mach 
number at shock formation, M = 0.01. 
The divergence may be remedied if it is recalled that the bulk 


absorptive processes have been neglected up to this point. If the 


viscous terms are retained and assumed to be additive to the wall losses 





then 
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The analysis and results are the same as before with the H's and 9's 


redefined as 


Be ie) G 
Hs - Cort Gu 4 





=e pt ida eee 
G, = Law Av be (4.5) 


With these modifications the amplitudes of the (n + 1)th order terms 


now behave, for sufficiently large harmonic, as 


PLT a = GS oF 7 


which tends to zero with increasing n. 

The modifications introduced by Eqs. 4.2 and 4.3 are important 
in the perturbation expansion for those harmonics whose indices 
Satisfy the inequality nef. . For the problem dealt with 
in this research the theory developed without the bulk absorptive 
terms is adequate as long as attention is restricted to frequencies 
below about the thousandth overtone. 

In comparing the theoretical and the experimental results it is 
convenient to divide the discussion into three areas: the low Mach 
number region, the intermediate Mach number region and the high Mach 
number region. 

In the low Mach number region, below M~0.006, it is obvious 
from Figs. C.1 and C.2 that theory and experiment are in very close 
agreement. In these figures of pressure waveforms the large drawing 
contains the waveforms predicted from the computer programs and the 
smaller inset is a drawing from an oscilloscope photograph of the 
waveform observed experimentally. The dashed line in both drawings 


is a sine wave which is included for reference purposes. One of the 
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main reasons for the good agreement in this region is due to the 
fact that the ("P/25. ) terms are as yet relatively small. 

It should be noted here that in all the figures of pressure 
waveforms, the major difference between the FINAMPI and FINAMPIZ 
predictions is in the phase dependence. Recall that the development 
of the reasoning behind the FINAMPIZ program rested on the possibility 
that corrective terms to the generated harmonics could be ignored. 
These curves lend support to that suggestion in the low Mach number 
region. 

The low Mach number region being discussed also includes the 
first two sets of Q-curves, Figs. C.5 through C.16. It is evident 
in both sets that the agreement of the experimental results with the 
predictions of the Fourier synthesis approach is excellent in the 
lower harmonics and deviates only slightly in the higher harmonics. 
When the agreement does deviate, the experimental results in general 
fall between the two theoretical predictions. This tends to indicate 
that the two solutions have the exact solution bracketed. It is 
evident in the QCURVES predictions for M = 0.005 that the third 
harmonic curve is no longer a smooth bell-shaped curve. The reason 
for this is the growth in the third harmonic of the correction terms 
arising from the fifth-order solution, For instance, the ratio of 
the approximate magnitude of the first corrective term to the leading 
term in the third harmonic is about 0.5 at M = 0.005. This is no 
longer insignificant so that the predictions based on the perturbation 


approach begin to show unrealistic fluctuations. 
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Theoretically predicted amplitude information is obtained from 
the two programs, PHAMP and PHAMPC. The amplitude results from these 
programs together with the experimental results are contained in 
Table 4,1. It should be noted that Beech's uncertainty of + 0.1 Hz 
in frequency results in an uncertainty of + 0.628 radians in Figs. (Gee 
through C.23. A comparison of the uncertainty in amplitude (based 
on the uncertainty in 4W) with the observed discrepancy between 


measurement and theory for M = 0.004 and4W= 0 yields 


n uncertainty in discrepancy 
harmonic content between PHAMPC 
and experiment 


Zz 0.5 0.4 

5 UES Orly 
4 OAD OFA 
5 Oo 0.08 
6 0.03 Or. 07 


It is evident that the uncertainties in harmonic content are consistent 
with the discrepancies in harmonic content noted between PHAMPC and 
experiment. 

As was stated in Section 3, PHAMPC also predicts the phase 
dependence of the various harmonics. The results of this program 
are contained in Table 4.2 and Fig. C.23. One thing should be noted. 
The sets of phase angles predicted by PHAMPC are identical for all 
Mach numbers, which is not surprising if the reader recalls that the 
normalization process destroyed the Mach-number dependence of the 


leading terms of each harmonic. 
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One other program exists for the investigation of phase and 
amplitude characteristics. The computer program FOUANAL is being 
introduced here since it was not derived from either theoretical 
approach. FOUANAL takes data from an oscilloscope photograph and 
computes the Fourier coefficients of the various terms of a Fourier 
coefficients of the various terms of a Fourier expansion of the 
waveform. The results from this program are also contained in 
Tables 4,1 and 4.7. 

It is clear that the results from FOUANAL are fairly accurate 
in the lower harmonics (first, second and third) but tend to become 
much worse in the higher harmonics (fourth, fifth and sixth). Judging 
from the amplitude results it does not seem likely that FOUANAL would 
yield accurate results in the phase angles of these higher harmonics. 
Examination of Table 4.2 verifies this supposition. These erratic 
results in the FOUANAL calculations may possibly be explained in the 
following manner. The input data for the program come from extracting 
sixty-four points from a standard oscilloscope photograph. This was 
done by hand and was subject to many errors. It is possible that the 
data were not given to sufficient accuracy to enable the computer to 
predict accurately past the third harmonic. The most feasible suggestion 
for improvement in this respect is to use a digital counting procedure 
to obtain the data. It should be noted that the data in Fig. C.23 
are the FOUANAL results for M = 0.004, Clearly these points are 
extremely random and the uncertainty in the experimental MW cannot 
possibly account for the entire discrepancy. 

Shifting now to the medium Mach number region (0.006 MS 0.009) 
we see from Fig. ©.3, for M = 0.007, that although the general agreement 


is still good, deviations from the previous good agreement are becoming 
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evident in the computer predictions. Unfortunately experimental 
data for Q-curve and phase angle comparisons were not obtained 

but it is not unreasonable to assume that these comparisons would 
show essentially the same result as stated for the pressure waveform. 

Moving into the region of high Mach numbers, M = 0.009, we find 
that the pressure waveforms predicted by each program are generally 
invalid. This is nearing the experimentally observed Mach number for 
shock formation (M = 0.01), and since the theory is predicted on small 
Mach number, pre-shock conditions it is not surprising to see these 
deviations. 

Essentially the same results are observed in the Q-curve comparisons, 
Figs. C.17 through €.22. The reason for the experimental amplitudes 
falling off so rapidly is at present unknown. Tables 4.1 and 4.2 
again reflect the serious discrepancies in the high Mach number region. 

It is now possible to make some general comments concerning these 
two approaches. Both the perturbation approach and Ene Fourier synthesis 
approach have adequately described the effects of finite-amplitude standing 
waves in rigid-walled cavities provided the Mach number is relatively low. 
This precludes that area in the pre-shock region where significant 
distortion exists. It has been demonstrated that in the region where 
the approaches are fairly accurate, Q-curve information can be easily 
obtained. This information is convenient for comparisons with experi- 
mental results. It is also possible to obtain some information on the 


phase angles of the harmonics. 
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It seems plausible now to suggest the following: Since the 
Fourier synthesis approach gave the most accurate information, the 
computer program involving all ten harmonics should be reviewed in 
the hope that a usable program will result. This program could be 
used to investigate the effects of the higher order harmonics on 
the existing solution and possibly give an indication of whether 
extending the present developments past six harmonics is beneficial. 

Finally, after some exposure to basic communications theory, it 
appears that a good grounding in this subject is mandatory if the 
present theoretical line is to be pursued further. It would not be 
surprising to find that a more exact solution may be obtained as a 


result of the techniques used in communications theory. 
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APPENDIX A 


Coppens-Sanders Iteration Scheme 





From Section 2, the approximate wave equation for finite-amplitude 


standing waves in a i walled tube is 


a ae Pr ron a 4 2) dng 
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For convenience the ro hand side of Eq. A.1 will be written 


2 a” 7 
Dx ——e : Age, cr (A. 2) 
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where the 7; and "7 are complex, and the three dots in the 

A/ 
exponential term represent phase terms to be defined. Work is done 
in one n at a time so that the n subscript may be dropped. 


Expansion of the left hand side of Eq. A.1 yields 
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where (-&) was previously defined to be Gr) 


Thus Eq. A.1 becomes 
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for a given order n. If Eq. A.5 is generalized in time it becomes 
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Subsestutionumer EquyA.3 imto"™Eq. A.6 yields 

Oz. 2? 2 «2 2 yp2 2 
= C7 byw =CG@ a k=) a ; G (AT) 
Cop k +s 7g D4 fy; 


Now each V will be of the form 
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Substitution yields 
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Additional substitution of Eq. A.9 and Eq. A.3 into Eq. A.7 gives 


UE) 4, -<f |- ¢ Anj fb “3 (A.10) 


Proper manipulation of the bracketed term in Eq. A.10 will yield the 


(A.9) 


phase relationships mentioned previously but left undefined. 


Fieste recall that 








a = a, WA 
' : (A.11) 
Op eh i of L sap i 
C Zs ™ 7 Ke 
where 
ha Me 
(Wa)? = t-G 
and notice that 
a yp tTOAW © x" Aw 
hn 2) eA (4.12) 
and 
SOs. 4 
> = Ge 2 4 (143 ae, (A,13) 


55 


ct alse gh + fa?) (A. 


Substitution into Eq. A.10 yields 


Gy (Goa BG) hy fg 


It is now possible to write down the now U as 
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Define the following quantities; 
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Pictorially these definitions may be represented as in Fig. A.1l. 





Spb + Lay), 


Pictorial representation of the quantities 
defined in Eqs. A.17 through A.19. 


Figure A.1 
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With these definitions, Eq. A.16 becomes 
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i Any GZ pe 





and 
Pr; “An Cox KE- a) cenljat ou) ve a (a2) 
yielding 


coy 4, 7 ng Cx) e t oep J (A.22) 


The preceding was a general outline of the procedure to be 
followed in the Coppens-Sanders iteration scheme. The remainder of 
this section contains a derivation of the terms through third order. 


Assume the input 


- bi Fi. Mea) Coa wt (A.23) 


which was shown in Section 2 to be the calculated first-order solution. 


Then we have 


£ = ve tin K{ka) hi OF (A. 24) 
OS, 
2 og De oop Hen) wwt a 


and since Oy - 
itinn. gh Ur 


a 





ee (LNW) a 
On Cs “7 <) Les at (A .26) 


The next step is to obtain 2° 


/ | 
f o) j n-¢ ' 


aa aa Oa das (A.27) 


Tae! 


B/ 


Substitution of Eq. A.26 into Eq. A.27 yields 


Wz =] Sane 2 MRL DHT. 


(A.28) 
where M = Un/C¢ 
Expansion of the trigonometric term gives 
Cor?*>Kb-a) ta *wTt (A.29) 


=A) fon feeteey} + [ow tul)-1] 


from which the constants + 1 and - 1 may be discarded since application 


2 
of oA dt will remove them. Then 


/ 
V2 = °F) Cor 2kC«) Poe Deol (A, 
which implies 
A. = / as 
2a al, (A. 
and from Eq. A.20 
We = ~My é mel 7 C CEz 
22 = ¢ A» pp . } 
- Ley De eo (A. 
2% (_ o ) pao 
Recall that Ang @ 4n7 , and hence, 
/ ME) is. 
42 = 5 Ulta AK L-a) Gor (20 +B2)n, 
which is the second order solution. 
The general equation for IZ a is 
) U/ A, 
om == 5 coe m Ma) shew (malt t+ 9), ( 
Substitution of Eq. A.34 into Eq. A.21 results in trigonometric 
such as 
Cr. mb(La) dient thy) Cos MK La) ti. lub + fy ) an 


Be 
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31) 


32) 


33) 


34) 


products 


35) 


For simplicity let 


eee 
T= va (A. 36) 


Then Eq. A.35 may be written 


Coa mX Li. MT Gy) Cor HX ta. (T+ te) (A. 37) 
Expansion of Eq. A.3/7 gives 


C— foo (mt 4)X Con (at) Tt by 4 G } 


(A. 38) 


g cet (nt) X Gnd ml) T + hd, f 


where terms of the form sin mK(L- a) cos Lot have been ignored 

in Eq. A.38 since substitution of these terms into the perturbed wave 
equation result in amplitudes of order tH which have been neglected 
throughout. Eq. A.38 must always be written as (m - 1) where m > i. 
since a positive frequency has always been assumed. 

Before proceding to the derivation of the third harmonic terms, 
a brief discussion might be of some use in acquainting the reader with 
higher order terms. 

In the third harmonic derivation, the user of this iteration 
scheme is first exposed to the development of more than one term in 
the solution for a particular order. In other words, the iteration 
yields both the leading term to the third harmonic plus the first 
correction term to the present solution. At the same time, the use 
of Eq. A.38 is demonstrated for the first time. It should also be 
pointed out that the methods used in the third harmonic iteration are 
identical to those used in the higher order iterations, regardless of 


the number of terms encountered. 
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Following the step by step procedure then, 


cf Cone FS AK OF? oi _ okt ahi 
Ws z, a an de ea Ja <n Ja 


Remember , 


- = os ti. Kl-a) bin. wl 


23 & —¢ be Cr2 Kb-a) hen wt 


i= ee (S4) x bi, 2A) Lieut + G2) 
oe ~2% LYE cor 2K E-4) yi lal Cb G2) 
= > (M. 7. ne: 
lie - £U,(B)4, a 


so that 


od fa Ls Coz KM L-4) hintl- dle 


cor ZR L-a) hw (ZC 4 6) } 


or @wanitig Eq. A.38, 


sg Ye Pj [few PR(L-a.) Caf Ful + G2) 


me Skea K(L-2) Cor (wt # ex) | 
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G39) 


(A.24) 


(A255 
(A .40) 
Cry 


(432) 


(A .42) 


(Aaa) 


Using Eq. A.8 implies 





Wi a ee Li Ue, 
PF vi Co Co (A.44) 





Ce CO; (A.45) 
and from Eq. A.10, 
- ) : 
Ue UR 2 fae (A.46) 
Zz Co 2 fis 


- UB) xe o § (G27 5) 
and 


U3, am gute C2 2 7 fae’? 
= -FUL Ma gae ((6,+ 62) 


Now Up = U2 se , and when the substitutions are made, the 


(A.47) 


final form of the third harmonic iteration becomes 


ether 


(A.48) 


= Z “il Bs Le KL-4) Cor(wt 4 6, +62). 
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A combination of the first three orders to obtain the total 


solution thus far results in 


a gym K(l-a.) Cn we 
"2s ($Y) jf, 44 2+) 001 Cat “ 
oth é v7 a) cea (Zut? G2) 


2 OW) iy, bi. HKl-s) 0m Cub 1b2 163) 
- $l) Ty, 4i. Ka) Cor (at + G4 62), 


Notice that the second term of the third order solution is of order 

one in j, (with reference to Eq. A.3), and thus it is termed the 

first correction term to the fundamental or first order solution. 
Rewriting Eq. A.49 in order to achieve a form consistent with 


this last observation and the form of Eq. 2.20 results in 


, 4 — eka) f Leo wt : 


—- (Mp) Colts i 46462 ) 
Fay f 


tik 2K(1-a) 13 (HF) J cor (att 62) 
a 


+i fie FMM SCE) Ly, Cor(Fat + Ge + @3) 
ist 


(A.50) 
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APPENDIX B 


Computer Programs 


An IBM-360 digital computer was utilized in performing the 
calculations, and obtaining the wave forms predicted by the theory. 
Computer programs were written for each of the basic formulations 
as well as for each of the variations investigated, Each of these 
programs, together with a short description, is presented in this 
section. Included with the description is a list of the symbology 
peculiar to that particular program. 

B-1l: Program FINAMPI 

This program has written to perform the calculations and plot 
the graphs for the perturbation solution. It was originally written 
in FORTRAN 60 language and subsequently converted when the IBM-360 
became available. 

The program was designed so that the Mach number, resonance 
frequency, length of the tube, position in the tube, 7, ba, and A, 
are the input parameters. The first half of the program is devoted 
to the velocity profiles and the second half to the pressure waveforms. 
The functions defined at the end of the main program are part of the 
conversion from FORTRAN 60 to FORTRAN IV. A list of the symbology 
pertinent to this program, together with the corresponding parameter 


in the original formulation follows. 
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Table B.l 


DIN) = & 
S = Mach number 
WD = Aw 
WR = Wy 
> =f 
F —ae 
A el 
U = Velocity 
Ie = Pressure 
PI = © 
aa ae 
.. 
Oi) Ss /2 Hp 
THET(N) = Gy 
Vxyz = the expansion of Oy), for the velocity equation 
VORA = the expansion of (X3y), for the pressure equation 
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B-2: Program FINAMPIZ 


This program was designed to duplicate the output of FINAMPI 
using the Fourier synthesis formulation through sixth order terms. 
Since the final results of the perturbation approach are in the same 
form as those of the Fourier synthesis approach the program FINAMPI 
was converted to FINAMPIZ by setting the corrective terms to zero. 


Because FINAMPIZ is essentially FINAMPI, Table B.1 is applicable here. 
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B-3: Program QCURVES 


This program was designed to compute and graph the normalized 


Q-curve for each of the six harmonics using the perturbation solution. 


In this program, the input parameters consist of the resonant 


frequency, the tube length, J P, and the piston acceleration. The 


symbology pertinent 
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F 
Delt(N) 
H(N) 
THET(N) 


SCRIPN(M) 


MAXA 


SMALA 


PHT 


to this program is given in the following table. 


il 


tl 


II 


Table B.2 


piston acceleration 

length of tube 

day 

hypothenuse of triangle described in Fig. A.1l 
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unnormalized magnitude of the n-th harmonic 
at point m 


the normalizing factor defined for each 
harmonic 
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B-4: Program QCURC 
This program is the Q-curve program based upon the Fourier synthesis 
approach using only the first six harmonics. The curves are normalized 


in exactly the same manner as in QCURVES and no new quantities have been 


defined. 
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_ he GRAPHS THE NORMALIZED Q CURVE 
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B-5: Program PHAMP 


This program was designed to compute the amplitudes and associated 
phase angles for each harmonic. Essentially this program is the QCURVES 
program with only minor modifications. One new quantity has been defined 
as PHASE(I,J). This is the phase angle associated with the i-th harmonic 


and the j-th Aa), In all other respects Table B.2 is applicable. 
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THF PROGRAM CUMPUTES THE AMPLITUODF AND PHASE OF 
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B-6: Program PHAMPC 


This program is identical to PHAMP except that it was written for 
the Fourier synthesis approach and graphs the phase angles associated 
with each harmonic. In order to use the graphing subroutine the 
subscripts on the phase angle had to be reversed. PHASE(I,J) is now 
defined as the phase angle associated with the i-th Aw and the j-th 


harmonic. In all other respects Table B.2 is applicable. 
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B-7: Program FOUANAL 


This program is an adaptation of work done by Roy M. Johnson, 
Assistant Professor, Naval Postgraduate School. Basically the program 
takes a waveform obtained directly from the cavity and performs a 
Fourier analysis on it. The output of interest is a set of Fourier 
coefficients and their associated phase angles. The program is 
written in such a manner that a table of definitions is not appropriate. 


Each variable is defined within the program. 
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Q-curves for the fundamental, M = 0.004, _ 
QCURC prediction,---QCURVES prediction 
Experimental results are indicated by @, 





Figure C.5 
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Q-curves-for the second harmonic, M = 0.004 
———QCURC prediction,----QCURVES prediction 
Experimental results are indicated by @ 


Figure C.6 
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Q-curves for the third harmonic, M = 0,.004,:- 


QCURC prediction,----QCURVES prediction 
Experimental results are indicated by @ 





Figure C,/ 
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Q-curves for the fourth harmonic, M = 0,004, 
QCURC prediction,----QCURVES prediction 
Experimental results are indicated by @- 





Figure C.8 
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Q-curves for the fifth harmonic, M = 0.004, 
QCURC and QCURVES prediction 
Experimental results are indicated by @. 





Figure os CC 
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Q-curves -for the Sixth harmonic, M = 0.004 
QCURC and QCURVES prediction 
Experimental results are indicated by@, 


Figure C.10 
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Q-curves for the fundamental, M = 0.005, . 
QCURC prediction,----QCURVES prediction 
Experimental results are indicated by @ 





Figure C.11. 
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Q-curves for the second harmonic; M = 0.005, 


QCURC prediction,----QCURVES prediction 
Experimental results are indicated by @. © 


Figure C.12. 
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Q-curves for the third harmonic, M = 0.005, 
QCURC prediction,---QCURVES prediction 
Experimental results are indicated by @- 
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Figure C.13 
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Q-curves for the fourth harmonic, M = 0.005, 
QCURC prediction,---QCURVES prediction 
Experimental results are indicated by @e 
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Figure C.14 
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Q-curve for the fifth harmonic, M = 0.005, | 
QCURC and QCURVES prediction 
Experimental results are indicated by @.. 
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Figure C.15 
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Q-curve for the sixth harmonic, M = 0.005, , 
QCURC and QCURVES prediction 
Experimental results are indicated by @o 


Figure C.16 
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Q-curves ‘for the fundamental, M = 0.009. 
- QCURC prediction,---QCURVES prediction 
Experimental results are indicated by@. 
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Figure C.17 
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Q-curves for the second harmonic, M = 0.009, 
QCURC prediction,----QCURVES ‘prediction 
Experimental results are indicated by @. 
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Figure C.18 
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Q-curves for the third harmonic, M = 0.009, 
QCURC prediction,----QCURVES prediction, 
Experimental results are indicated by @, 


Figure C.19 
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Q-curves for the fourth harmonic, M * 0.009, 
QCURC prediction,---QCURVES prediction 
Experimental results are indicated by @, 


Figure C.20 
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Q-curve for the fifth harmonic, M = 0.009, 
QCURC and QCURVES prediction 
Experimental results are indicated by @. 
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Q-curve for the sixth harmonic, M = 0.009, 
QCURC and QCURVES prediction 
Experimental results are indicated by, 


Figure C.22 
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Phase dependence of the various harmonics on A@® 
as predicted by PHAMPC. n = harmonic number, @ 
= data from FOUANAL for n'-th harmonic 


Fi, * 93 


145 


INITIAL DISTRIBUTION 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library 
Naval Postgraduate School 
Monterey, California 93940 


Commander 


LIST 


No. Copies 


20 


Naval Ships Systems Command Headquarters 


Washington, D. C. 20360 


Commander 


Naval Ordnance Systems Command Headquarters 


Washington, D. C. 20360 


Professor Alan B. Coppens 
Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 


Professor James V. Sanders 
Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 


CDR Wayne "'L'" Beech 
1108 Leahy Road 
Monterey, California 93940 


LY al Geakurt it 


USS Dale (DLG~19) 
FPO San Francisco, California 


146 





UNCLASSIFIED 
Security Classification 


DOCUMENT CONTROL DATA - R&D 


(Security classification of title, body of abstract and indexing annotation must be entered when the overall report is classified) 


Mie! Fosteraduate School 
Une lasstit tea 


Naval Postgraduate School 
Monterey, California 93940 Sb anaoe 









3. REPORT TITLE 


A Theoretical Investigation of Finite Amplitude Standing Waves in Rigid Walled 
Cavities 





4. DESCRIPTIVE NOTES (Type of report and inclusive dates) 


Master's Thesis, December 1967 
5. AUTHOR(S) (Last name, firat name, initial) 


mePe, Paul G., IIL 


6. REPORT DATE Ja. TOTAL NO. OF PAGES 7b. NO. OF REFS 


December 1967 | 144 =26 


8a. CONTRACT OR GRANT NO. 9a. ORIGINATOR’'S REPORT NUMBER(S) 


b. PROJECT NO. 


y 10. AVAILABILITY/LIMITATION NOTICES 


a ae eon | = . | 


Oo 


11. SUPPL EMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
Monterey, California 93940 


13. ABSTRACT 


The Coppens-Sanders perturbation solution for the one-dimensional non- 
linear acoustic wave equation with dissipative term describing the viscous 
and thermal energy losses encountered in a rigid walled, closed tube with 
large length-to-diameter ratio was extended to include sixth order terms. 

The solution was then investigated to determine the region of validity. 
Computer programs were written to evaluate and graph the resulting waveforms. 
Available experimental results were compared with the theoretical predictions 
and good correlation was found to exist in the region of low Mach numbers. 
This agreement was found to gradually deteriorate as the Mach number was 
increased. A Fourier synthesis approach is also presented and the leading 
terms of the first ten harmonics are derived. 





DD aon 6 1473 147 UNCLASSIFIED 


Security Classification 


UNCLASSIFIED 


Security Classification 


14. 
aye LE Nie A LINK #8 LINK ¢ 





Finite-AmpLlitude 

Standing Waves 
Perturbation Approach 
Fourier Synthesis Approach 
Rigid-Walled Cavity 


% 
¢ iia - oe 


ae 3 $ = 3 
we 


FORM | 
DD ev 1473 (Back) 148 UNCLASSIFIED 


S/N 0101- . Ge fall PG en: Sen an ae 
1-807-6821 Security Classification A-31409 











‘main 


-_ 2768 00414127 5 





